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We investigate electronic transport through parallel double quantum dot(DQD) system with 
strong on-site Coulomb interaction and capacitive interdot coupling. By applying numerical renor- 
malization group(NRG) method, the ground state of the system and the transmission probability at 
zero temperature have been obtained. For a system of quantum dots with degenerate energy levels 
and small interdot tunnel coupling, the spin correlations between the DQDs is ferromagnetic, and 
the ground state of the system is a spin 1 triplet state. The linear conductance will reach the unitary 
limit (2e 2 /h) due to the underscreened Kondo effect at low temperature. As the interdot tunnel 
coupling increases, there is a quantum phase transition from ferromagnetic to anti-ferromagnetic 
spin correlation in DQDs and the linear conductance is strongly suppressed. 

PACS numbers: 72.15.Qm, 73.23.Hk, 73.40.Gk 



I. INTRODUCTION 

In recent years considerable research attention has 
been paid to electron transport through double quan- 
tum dot (DQD) systems p], which are an artificial small 
quantum systems can be readily controlled by exter- 
nal gate voltage and also exhibit a variety of inter- 
esting strongly correlated electron behaviors. Basi- 
cally, there are two different experimental realizations of 
DQD system: DQDs connected in serial[2J or in parallel 
configurations [5] ■ Electron transport through both con- 
figurations have been studied in experiments, the molec- 
ular states of the double dots and also the competition 
between Kondo effect and the RKKY interaction have 
been observed [2J [3]. 

The theoretical studies on electron transport through 
DQDs are largely devoted to the system in the Kondo 
regime. For DQDs connected in serial, the antiferro- 
magnetic correlations between two single-level coupled 
QDs is in competition with Kondo correlations between 
the QDs and the electrons in the leads, therefore it 
gives rise to rich ground state physical properties at 
zero temperature [U [5J |5J [7]. For DQDs with large ca- 
pacitive coupling, the simultaneous appearance of the 
Kondo effect in the spin and charge sectors results in 
an SU(4) Fermi liquid ground state[S]. By increasing in- 
terdot capacitive coupling, a quantum phase transition of 
Kosterlitz-Thouless type to a non-Fermi-liquid state with 
anomalous transport properties is predicted[9]. Martins 
et al argued that ferromagnetic state cannot be realized 
in two single-level QDs connected in serial, but they pre- 
dicts that FM state can be developed in two double-level 
QDs [10 . For the DQD system in parallel configuration, 
the physical properties can be quite different, since the 
interference effect will play an important role in its trans- 
port properties. The Fano effect for electron transport 
through bonding and antibonding channels in DQDs sys- 
tem has been studied [TT1 H"2l IT5] . 



Due to the strong correlation of electrons in the QDs, 
it is a non-trivial problem to treat those systems theo- 
retically. It is well known that Wilson's numerical renor- 
malization group [121 HS1 HSJ EZ] method is a nonpertur- 
bative approach to quantum impurity problem, which 
can take into account the on-site Coulomb repulsion and 
the spin exchange interaction between the electrons in 
DQDs exactly, in contrast to the slave boson mean field 
theory or the equation of motion method within Hartree- 
Fock approximation. The NRG method have already 
been applied to investigate a lot of problems in the elec- 
tron transport through QD systems. For instance, DQDs 
connected in serial [15]. the quantum phase transition in 
multilevel QD[19 , Kondo effect in coupled DQDs with 
RKKY interaction in external magnetic field [20 , the side 
coupled DQD system[2f| [55] and quantum phase transi- 
tions in parallel quantum QDs |23j etc. However, in our 
opinion the consequences of the interplay of Fano res- 
onance and the Kondo effect on electron conductance 
through DQDs in parallel still haven't been well eluci- 
dated. In this paper we will investigate the electron 
transport properties for the DQDs in parallel configu- 
ration by using the NRG method. We will show that 
for DQDs without interdot tunneling, the underscreen 
Kondo effect plays an essential role in the conductance. 
The linear conductance, spin correlation, and local den- 
sity of state in this system are obtained. 



II. THE MODEL HAMILTONIAN AND THE 
NRG APPROACH 

Electron transport through parallel-coupled DQDs 
with interdot tunneling, on-site Coulomb interaction and 
capacitive interdot coupling can be described by the fol- 
lowing Anderson impurity model: 
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where cjfe^ ff (c^ ) denote annihilation (creation) operators 

for electrons in the leads(?7 = L,R), and those 
of the single level state in the i-th dot( i = 1,2). m a 
denotes the electron number operator with spin index a 
in the i-th dot, and rii — ^2 a ni a . U is the intra-dot 
Coulomb interaction between electrons, V is the interdot 
capacitive coupling. t c is the interdot tunnel coupling, 
and v v i is the tunnel matrix clement between lead r\ and 
dot i. It should be noted that an interdot magnetic ex- 
change term J is not explicitly included in this Hamil- 
tonian since it is not an independent parameter but a 
function of the interdot tunneling (J ~ t^/U). We con- 
sider the symmetric coupling case with Tf = Tf = T{ . 
where = 27r^ fe \Vm\ 5{u> — e^a) is the hybridization 
strength between the i-th dot and the lead 77. 

In order to access the low-energy physics of this DQD 
system, we adopt the Wilson's NRG approach. By sym- 
metric combination of the lead orbitals, the Hamiltonian 
in eq.(l) can be mapped to single-channel two-impurity 
Anderson model. Because the anti-symmetric combina- 
tion of lead orbitals are totally decoupled with the QDs, 
they can be neglected in the Hamiltonian. Following the 
standard NRG method, one defines a series of rescaled 
Hamiltonian Hn as following 
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Green's functions(GF) G r /°(w) have 2x2 ma- 
trix structures, which account for the double dot 
structure of the system. The matrix elements 
of the retarded GF are defined in time space as 
G^t-t 1 ) = -i6(t-t') < {d 4(T (i),d+ (i')} >. Therefore, 
the transmission probability T(u>) can be obtained by 
calculating the imaginary parts of the GF of DQDs or 
the spectral density pij{u) = — -7mGy(u;). Then, the 
linear conductance at the absolute zero temperature 
can be given by taking the zero frequency limit of the 
transmission probability G = J^|y = o = ^f^T(u) = 0). 
One advantage of the NRG is accurate determination of 
the low-energy spectral density of the quantum impurity 
models. By a standard procedure in NRG [17], the 
spectral density at zero temperature can be calculated 
according to the following formula 



PiM = ^^Ml x (Ml x yS(uj-(E x -E Q )) 

+ ^5yE M A.o( M lo)*^ + (^A-i?o))(5) 



where the matrix element M\ =< A|d,v|0 >, with |0 > 
and I A > being the ground state and exited eigenstate of 
the impurity model Hamiltonian, respectively. 
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where the discretization parameter A 
l|17j. The other parameters e, = 
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( 1+A.- 1 ) 2 ^b ' w ^ n D being the bandwidth of electrons 
in the leads. The above one dimensional lattice model is 
iteratively diagonalized by using the recursion relation 

Hn+1 — A}/ 2 H N + £ n ^^{Sl] a fN+la + flq + i a fNa) ■ (3) 



The basis set in each iteration step is truncated by retain- 
ing only those states with low-lying energies. In our nu- 
merical calculation, we keep totally 600 low-lying energy 
states in each step without counting the S z degeneracy. 

The current formula through the DQDs is given by the 
generalized Landauer formula [24 



(4) 
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FIG. 1: (a) The electron occupation number < m > in each 
quantum dot as a function of the gate voltage e<j- Ae^/r = 
0(solid line); 2.0 (dashed line). The other used parameters 
are D = 1.0, t c = 0, F = 0.01, U/T = 10, and V = U/2; (b) 
The interdot spin correlation < Si • S2 > as a function of the 
gate voltage e^. 
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III. RESULTS AND DISCUSSIONS 

In the following, we will present the results of our 
NRG calculation. For the sake of simplicity, we only con- 
sider the symmetric coupling case with the hybridization 
strength r a = T 2 = T. We take the bandwidth D = 1 
as the energy unit, and the other parameters T = 0.01, 
U = 10T, V — U/2. One can define the averaged energy 
level of QDs as = (ei + £2)/2, and the energy level 
difference Ae^ = e 2 — t\. Both of them can be tuned 
experimentally by external gate voltages. 

At first, we consider DQDs without interdot tunneling 
(t c = 0). In Fig. 1(a) the occupation number of electrons 
< rii > in each QD is plotted as a function of the av- 
erage energy level e^. The electron occupation number 
increases consecutively by tuning the QD level below the 
Fermi energy. For this DQDs with interdot capacitive 
interaction, one can easily discern the different regions of 
occupation states: from empty occupation to the state 
with total four electrons in DQDs . In the case of two 
identical QDs ( Ae^ = 0), abrupt jumps of the occupa- 
tion number are observed at some particular gate voltage. 
One can see that the position of jumps can be identified 
as the region where the DQDS have odd number of elec- 
trons. For DQDs with different energy levels(Aed ^ 0), 
the QD with low energy level is occupied first, and be- 
cause of interdot capacitive interaction, it will greatly 
suppress the occupation of electron in another QD as 
compared with the two identical QDs case. The inter- 
dot spin-correlation < Si • S2 > as a function of energy 
level Cd is shown in Fig. 1(b), where the spin operators in 
the i-th QD are defined by Si = 1/2 d\ a a acr 'di a r . It 
shows that the interdot spin correlation is antiferromag- 
netic in the mixed valence regime, and is ferromagnetic in 
the doubly occupied regime, where each QD is occupied 
by one electron. For DQDs with energy level difference, 
the spin correlation in mixed valence regime is greatly 
suppressed, but there are still large ferromagnetic spin 
correlation in the doubly occupied regime. In the identi- 
cal QDs(Ae^ = ) case, the abrupt jumps in occupancy 
and the spin correlation turn from FM to AFM have also 
been found in Ref.gS] for iV-QD system(A > 2) with- 
out interdot capacitive coupling, and this phenomenon 
is interpreted as a kind of quantum phase transitions. 
However, one can see from the dashed line in Fig.l that 
the abrupt jumps both in occupancy and spin correla- 
tion disappear when Ae^ ^ 0, hence this kind of phase 
transition is unstable with respect to the perturbation by 
gate voltage difference in QDs. We attribute this kind of 
abrupt jump as a result of Fano resonance and the cross- 
ing of the antibonding state energy level with the Fermi 
energy. 

Next, we calculate the electron conductance through 
DQDs when a small bias voltage is applied to the leads. 
In Fig. 2 the linear conductance G at zero temperature 
vs. the average QD energy level is depicted. As shown 
in Fig. 2(a), the Kondo effects are manifested by peaks in 
the curve of the linear conductance, where the conduc- 
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FIG. 2: The linear conductance G as a function of the dot 
level at zero temperature, (a) for the system with two identi- 
cal quantum dots(Ae<j = 0); (b) for DQDs with energy level 
difference(Ae,j/r = 2.0). The interdot tunneling parameter 
takes t c /r = 0(solid line) and t c /T = 2.0(dashed line), re- 
spectively. 

tances approach the unitary limit (G — 2e 2 /h). In the 
regime of odd electron occupation, the DQDs act as a lo- 
calized spin (s=l/2), and the Kondo effect is arose from 
the spin exchange interaction between the localized elec- 
tron spin and that of the electrons in the leads. Whereas, 
in the doubly occupied regime, the unitary conductance 
is due to the underscreened spin 1 Kondo effect, It will 
be shown in the following that in this regime the two 
electrons confined in the DQDS form a spin triplet state 
in the ground state. In the presence of sufficient inter- 
dot tunneling t c , the Kondo effect in the singly occupied 
regime and spin 1 Kondo effect is strongly suppressed, 
but some asymmetrical peaks of conductance appears in 
the mixed valence regime, this can be attributed to the 
Fano resonance for the electron transport through the 
bonding and antibonding channels in this system. It is 
interesting to notice that in the triply occupied regime 
the conductance still achieves the unitary limit even in 
the presence of strong interdot tunnel coupling. 

In the following, we will focus our attentions on the 
properties in the doubly occupied regime. In order to il- 
lustrate the effect of interdot tunneling, the transmission 
probability at different tunneling coupling t c is shown in 
Fig. 3 (a). Without direct interdot tunneling (t c = 0), 
one can see that the transmission probability has the 
particle-hole symmetry, and the spin exchange effect be- 
tween the electrons localized in the quantum dots and 
that in the leads gives rise to a sharp peak in the trans- 
mission probability at the Fermi surface, therefore the 
linear conductance at zero temperature reaches the uni- 
tary limit G — 2e 2 /h, as a result of the underscreened 
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FIG. 3: (a) The transmission probability T(lo) for the sys- 
tem with two identical quantum dots. Used parameters are 
D = 1.0, T = 0.01, U/F = W,V= U/2 and e d /F = -10.0. 
The interdot tunnel coupling t c /Y = 0.0,0.8, 1.0, 1.5, respec- 
tively, (b) The transmission probability T(ui) at the particle- 
hole symmetric point for different values of on-site Coulomb 
interaction U. Inset: the estimated Kondo temperature Tk 
vs. the value of U. 



spin 1 Kondo effect. In the presence of the interdot cou- 
pling t c ^ , the particle-hole symmetry of the trans- 
mission probability is broken. When t c increases beyond 
a quantum critical point, a sharp dip in the transmission 
probability is observed. It suggests that the Kondo ef- 
fect and the linear conductance in this regime is strongly 
suppressed. Therefore, there is a quantum phase transi- 
tion between underscreened Kondo phase and the local 
spin singlet phase in the ground state of this system. For 
two-impurity Anderson model without interdot capaci- 
tive coupling, this quantum phase transition has been 
predicted by Nishimoto et al. by using dynamic den- 
sity matrix renormalization group [25], and Zitko et al. 
have obtained its thermodynamic properties, such as the 
temperature dependence of magnetic susceptibility and 
entropy by NRG method It is noted that a similar 
quantum phase transition is also observed in two-level 
single QD system with intradot spin exchange coupling 
by Hund's rule 19J. For DQDs with RKKY interaction 
coupled to two-channel lead, Chung et al-J^O] found the 
quantum phase transition is from Kondo screened phase 
to spin singlet phase. In Fig. 3 (a), by further increasing 
the interdot coupling t c , a broad peak of transmission 
probability with the line shape of Brcit-Wigner resonance 
is developed around the energy lu U/2, we attribute 
this broad transmission peak to the electron transport 
through the bonding channel of electrons in the quan- 
tum dots. 

In Fig. 3(b) the transmission probability T(ui) at differ- 



ent values of on-site Coulomb interaction U is depicted. 
It shows that the line shape of the T(uS) changes signif- 
icantly by varying the Coulomb interaction strength U . 
The line shape becomes more cusplike with decreasing 
U , and it reveals that the physical properties of this 
underscreened Kondo effect in DQD system is quite dif- 
ferent from the spin 1/2 Kondo effect. For the spin 1/2 
Kondo effect in single- impurity Anderson model, one can 
estimate the Kondo temperature Tk by using the for- 
mula T K = ^^exp[ed(e d + U)/UT]. For this under- 
screened Kondo effect case, we make the following ap- 
proximation to estimate the Kondo temperature : at 
the frequency of u> = Tk the transmission probability 
T(lu = T k )/T{uj = 0) fa 0.978. For the single-impurity 
Anderson model, Tk obtained by this approximation 
agrees well with the above formula. The inset of Fig. 3 (b) 
shows the estimated Tk at several values of the Coulomb 
interaction strength U for the DQD system. For the 
system with the parameters used in our calculation, the 
Kondo temperature Tk is on the order of 10 _5 r. 
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FIG. 4: The density of state of the local bonding and an- 
tibonding states in the quantum dot at different values of 
energy level difference: Aed/F = 0.0(solid line), 2.0(dashed 
line), 4.0(dotted line), (a), (b) corresponds to the bonding 
state with t c /r = 0,0,2.0, respectively. (c),(d) are that of 
the antibonding states. The other used parameters are the 
same as in Fig. 3. 

In order to get better understanding of the electron 
state in the system, we investigate the local density of 
states(DOS) in the DQDs. One can define the even or- 
bital (bonding state) operator as d ejCr = (di a +d 2 o-)/v^ ; 
and the odd orbital (antibonding state) operator rf e cr = 
(die — d2a-)/v2- The local density of state for the bond- 
ing and antibonding states are depicted in Fig. 4. As 
shown in Fig. 4(a) and (c), in the absence of interdot 
coupling(t c = 0), the local DOS of even and odd orbital 
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retain the particle-hole symmetry of the system. It is no- 
ticed that the transmission probability is proportional to 
the DOS for the bonding state, therefore a Kondo peak 
around the Fermi energy is observed in its DOS. Some 
new features are also manifested in DOS for this system, 
one can see that the DOS for the antibonding state has 
two side peaks nearby the Fermi energy , which can be 
understood as a result of the effective spin-exchange in- 
teraction between the electrons in DQDs by tunneling 
through the leads, and this feature cannot be found in 
DQDs in serial configuration [T5]. As the interdot cou- 
pling t c is large than some critical value(see Fig. 4(b) 
and (d)), the Kondo effect on the DOS of bonding state 
is greatly suppressed, and a broad peak around the en- 
ergy u> ss U/2 are developed. For the DOS of the anti- 
bonding state, a sharp peak is developed slight below the 
Fermi energy, this is due to the fact that the antibonding 
state of electrons in DQDs seems like a quasi-localized 
state. Increasing the interdot coupling t c further, the 
sharp peak is broaden and shifts away from the Fermi 
surface to lower energy. For DQDs with different energy 
levels, the characteristic features of the DOS remain un- 
changed. 
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FIG. 5: (a)The spin correlation < Si ■ S2 > of double quan- 
tum dots as a function of temperature T for several different 
values of t c . (b) The spin correlation < Si ■ S2 > vs. the 
interdot tunnel coupling t c at zero temperature. The other 
used parameters are the same as in Fig. 3. 

To gain more insight into the spin entanglement and 
the effect of spin exchange interaction for the electrons 
localized in different QDs, We have also calculated the 
interdot spin-correlation < Si • S2 >asa function of 
temperature for several values of interdot tunneling t c as 
shown in Fig. 5(a). When interdot tunneling t c is zero or 
has a small value, one can see that the spin-correlation 
converges to a positive value as temperature decreases. 
It is easy to notice that the positive value of < Si • S2 > 



reveals that the spin-correlation in this case is ferromag- 
netic type in the ground state. As we know that when 
two ideal spin s — 1/2 electrons form a spin triplet, the 
spin-correlation will be < Si-S2 >= 1/4. The rather high 
positive value of spin-correlation indicates that electrons 
localized in QDs still have high probability to form a 
spin triplet even though they are coupled with the elec- 
trons in the leads in the Kondo regime. By increasing 
the interdot coupling t c , there exhibits a quantum phase 
transition from the triplet state to singlet state in the 
ground state. The spin correlation approach a negative 
value < Si • S2 >~ —0.50, as we know that for two elec- 
trons forming an ideal spin singlet< Si • S2 >= —0.75, 
therefore the electrons in DQDs is largely in a singlet 
state. In order to determine the critical value of t c , we 
have calculated the spin correlation < Si • S2 > at zero 
temperature for different value of t c , the result is shown 
in Fig. 5(b). We find that, at the quantum critical point 
t c ra 0.7, there is an abrupt jump of the spin correlation 
< Si • S2 >• It indicates that the quantum phase tran- 
sition from triplet to singlet state is of first order kind. 
According to a previous study on two-impurity Kondo 
model |2"7] . we may expect that in the case of DQDs with 
energy level difference, this kind of first order transition 
will become Kosterlitz-Thouless type. It is easy to under- 
stand that the exact quantum critical value of t c will de- 
pends on the the interaction parameters, such as the on- 
site Coulomb repulsion U, interdot capacitive coupling V 
and the energy level ea, etc. 



IV. SUMMARY 



In summary, we have studied the ground state and the 
electron transport properties of the system with DQDs 
in parallel configuration. The strong on-site Coulomb 
repulsion and the interdot capacitive coupling is taken 
into account by the nonperturbative NRG technique. It 
is shown that the large interdot tunneling will drasti- 
cally change the transport properties in this system. The 
ground state of DQDs exhibits a quantum phase transi- 
tion from triplet state to singlet state by increasing the 
interdot tunneling amplitude. In the case of no interdot 
tunneling, the linear conductance approaches to unitary 
limit in the doubly occupied regime due to the under- 
screened Kondo effect, whereas it is greatly suppressed 
when the electrons in DQDs form an singlet state with 
the interdot coupling t c being larger than the critical 
value. For the DQDs with strong interdot tunneling, the 
Fano resonance can be observed in the linear conductance 
when the system are in the mixed valence regime. One 
may expect that the underscreened Kondo effect can be 
observed in future experiments on DQD system without 
direct interdot tunneling. 
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